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>r~i . Abstract 

We investigate the Casimir force between a microfabricated elliptic cylinder (cylindrical lens) 

o\ 

and a plate made of real materials. After a brief discussion of the fabrication procedure, which 
typically results in elliptic rather than circular cylinders, the Lifshitz-type formulas for the Casimir 
force and for its gradient are derived. In the specific case of equal semiaxes, the resulting formulas 

X 

coincide with those derived previously for circular cylinders. The nanofabrication procedure may 
also result in asymmetric cylindrical lenses obtained from parts of two different cylinders, or rotated 
through some angle about the axis of the cylinder. In these cases the Lifshitz-type formulas for the 
Casimir force between a lens and a plate and for its gradient are also derived, and the influence of 
lens asymmetry is determined. Additionally, we obtain an expression for the shift of the natural 
frequency of a micromachined oscillator with an attached elliptic cylindrical lens interacting with 
a plate via the Casimir force in a nonlinear regime. 

PACS numbers: 31.30.jh, 12.20.Ds, 12.20.Fv, 77.22. Ch 
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I. INTRODUCTION 



In recent years the Casimir effect 1] is acknowledged to be among the most rapidly 
developing fields of fundamental physics. It has attracted considerable attention as a test 
for the structure of the quantum vacuum, and for hypothetical interactions predicted in many 
extensions of the standard model, and also opened up new opportunities for nanotechnology 



12] 



Since 1997 approximately 30 experiments on measuring the Casimir force have been 
performed (see reviews 0, 4]), which not only confirmed the currently available theoretical 
knowledge, but also led to unexpected results of major importance. Specifically, it was 
recognized j^J that the unified theory of van der Waals and Casimir forces developed by 
Lifshitz encounters problems in the description of free charge carriers. As a result, two 



theoretical approaches were proposed based on the Drude 



-[8| and plasma |9Mll| models. 



Lifshitz theory, combined with the seemingly most natural Drude model, was shown to be 



in contradiction with the Nernst heat theorem 0, [12I 14] and with the experimental data 
In contrast Lifshitz theory using the plasma model for the dielectric permittivity 
was found to be thermodynamically and experimentally consistent, despite the fact that it 
does not take into account the relaxation properties of free charge carriers. Note that the 
experiment of Refs. [3, is an independent measurement of the gradient of the Casimir 
force with no fit ting parameters, such as a distance offset, etc. In the first repetition of 
this experiment [17| measurements were performed at separations up to 1.15 /im where 



zero values of the force were achieved within the limits of experimental errors. It was 



shown 



18| that the introduction of an offset did not improve the agreement between data 



and the Drude model and made the agreement with the plasma model worse. Because 
anomalous electrostatic contributions were not observed, introducing additional parameters 
in the theoretical description of the experimental data was unwarranted. 

This situation has been the subject of much controversy (see, for instance, Refs. 19h23|). 
Along with the experimental and theoretical investigations mentioned above, great progress 
was achieved in the calculation of the Casimir force between nonplanar surfaces based on the 



scattering approach [2I. I24L 29j . Bearing in mind, however, that in the end the elements of 
a scattering matrix are expressed in terms of dielectric permittivity or some other quantity 
characterizing material properties of the test bodies, successful application of new methods 
calls for the resolution of the problem of free charge carriers. 
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Presently great interest is expressed in new measurements of the Casimir force which 



could shed light on this problem. Thus, the experiment j30j claims observation of the 
thermal Casimir force, as predicted by the Drude model approach, in the separation region 
from 0.7 to 7/zm. It should be mentioned, however, that in Ref. 30] what is measured is 
not only the thermal Casimir force, but up to an order of magnitude greater total force 
presumably determined by large surface patches. The theoretical expression for the total 
force contains two fitting parameters determined from the best fit between the experimental 



data and theory. Therefore, Ref. 
force as is the experiment of Refs. 



30| is not an independent measurement of the Casimir 



In addition, it was shown 



3l| that the simplest 



version of the proximity force approximation (PEA), used in Ref. 30[ to calculate both the 
Casimir and electric force between a spherical lens with R = 15.6 cm radius of curvature and 
a plate, is inapplicable to large lenses due to the presence of surface imperfections. Another 
recent experiment employing large spherical lenses 32) does not support the existence of 
a large thermal correction to the Casimir force predicted by the Drude model approach. 
Because of this, new experiments, especially exploiting more sophisticated configurations 
than a sphere or a spherical lens above a plate, may lead to more reliable results than those 



obtained in Refs. 



30 



32]. 



As a prospective alternative configuration for the measurement o: 



33 



the Casimir force, a 



-4351] . This geometry is 



cylinder-plate geometry has long been discussed in the literature 
intermediate between the configurations of two parallel plates and a sphere above a plate. 
It preserves some advantages of the latter while making the problem of preserving the 
parallelism less difficult than for two plates. However, the configuration of cylinders with 
centimeter-size radii of curvature revealed anomalies in electrostatic calibrations 35] which 



might be caused by surface imperfections. To avoid this problem, Ref. [36| proposed an 
experiment measuring the thermal Casimir interaction between a plate and a microfabricated 
cylindrical lens attached to a micromachined oscillator. Such metallic lenses, with smooth 
surfaces of about 100 /mi radii of curvature on top of a micromachined oscillator, can be 
directly fabricated by using a monolithic fabrication process. In Ref. 36[ the Lifshitz-type 
formulas for the thermal Casimir force between a circular cylinder and a plate made of real 
metals were derived using the PEA. From a comparison with exact results available for ideal 
metals it was shown that for reasonable experimental parameters the error resulting from the 
use of the PFA is much less than 1%. This conclusion was confirmed in Ref. 37J for an ideal 
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metal cylinder above an ideal metal plate. It was shown that in the region of experimental 
temperatures the PEA correctly reproduces the dominant contributions to both the Casimir 
force and thermal correction to it. The validity of the PFA was also confirmed [38| for the 
configuration of an atom near an ideal metal cylinder. Reference demonstrated the 
feasibility of the proposed experiment, and investigated corrections to the Casimir force and 
its gradient due nonparallelity of a plate and a cylinder axis, and due to the finite length of 
a cylinder. 

In this paper we investigate the Casimir force between a microfabricated elliptic cylinder 
and a plate using the PFA approach. Our consideration is adapted to the measurement 
scheme using a micromachined oscillator. Motivation for use of elliptic cylinders derives 
from the fact that fabrication procedures usually result in cylinders with semiaxes in two 
perpendicular directions varying by 20%-30%. Fabrication may result also in asymmetric 
cylindrical lenses consisting of parts of two different elliptic cylinders or rotated through 
some angle about the cylinder axis. In all these cases we derive the Lifshitz-type formulas 
for the Casimir force and for its gradient and perform computations to account for the role 
of asymmetry. The electric force between an elliptic cylindrical lens and a plate is also 
calculated for the purpose of electrostatic calibration of the Casimir setup. Furthermore, we 
consider an elliptic cylindrical lens attached to a micromachined oscillator and interacting 
with a plate via the Casimir force in the dynamic regime. We derive the exact expression 
for a shift of the natural frequency of the oscillator under the influence of the Casimir force. 
This allows measurements of the frequency shift in a nonlinear regime, and comparison of 
the experimental results with different theoretical approaches to the Casimir force. 

The paper is organized as follows. In Sec. II the experimental procedures for microfab- 
rication of smooth cylindrical objects are considered. In Sec. Ill the Lifshitz-type formulas 
are derived for the Casimir force and for the gradient of the Casimir force between a plate 
and an elliptic cylinder. In Sec. IV the same is done for an asymmetric cylindrical lens near 
a plate. Section V is devoted to the micromachined oscillator with an attached cylindrical 
lens under the influence of the Casimir force in a nonlinear regime. Section VI contains our 
conclusions and discussion. 
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II. TECHNIQUES FOR MICRO FABRICATION OF SMOOTH CYLINDRICAL 
OBJECTS 



There are several approaches to creating a smooth cylindrical object on top of a micro- 
machined torsional oscillator that are fully compatible with ion chromatography techniques. 
As a consequence, the cylindrical microstructures can be monolithically integrated with the 
microelectromechanical oscillators. Examples include Focused Ion Beam technology (FIB) 
and techniques based on femtosecond laser microfabrication. 

FIB technology uses a Ga + ion beam to remove material from almost any surface. The 
profile to be patterned can be automatically inputted and controlled rather precisely. These 
tools are available in almost any microfabrication laboratory and, when combined with 
a scanning electron microscope, they can be used for non-destructive imaging at higher 



magnifications, permitting extremely accurate control of the milling process [39|,|40|. Today's 
most advanced FIB tools allow direct patterning of metals with minimum contamination or 
damage, which opens up the possibility to directly pattern nanostructures with desirable 
shapes onto microelectromechanical oscillators. 

Techniques based on femtosecond laser microfabrication, similar to the ones used to 



fabricate microlenses on glass 



4l|, can be also used to integrate a cylinder onto the paddle 



of a microelectromechanical system. In this case a tightly focused femtosecond laser beam 
is scanned inside a photosensitive material to create the required shape as precisely as 
possible. Once the photosensitive material is developed, the exposed volume will remain on 
the oscillator plate and standard etching processes can be used to transfer this shape to the 
plate. 

These microfabrication techniques represent, in our opinion, the most versatile methods 
of fabricating microstructures of a desirable shape. For objects of cylindrical shape, micro- 
fabrication typically results in elliptic rather than circular cylinders. The actual shape of a 
microfabricated object can be measured very precisely using a noncontact optical profilome- 
ter. Microfabricated cylinders may have semiaxes in two perpendicular directions varying 
by 20%-30%. They might be also characterized by some asymmetry (for instance, the axis 
of a cylinder may be not exactly parallel to the plate of a microelectromechanical oscilla- 
tor). More complex and expensive fabrication methods can also be used depending on the 



precision, uniformity and reproducibility required 



42 



43|. 
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III. THE CASIMIR FORCE BETWEEN AN ELLIPTIC CYLINDER AND A 
PLATE WITHIN THE PROXIMITY FORCE APPROXIMATION 



We consider an elliptic cylindrical lens of thickness h and width 2d obtained from an 
elliptic cylinder made of a material with a frequency- dependent dielectric permittivity e(u). 
Let the surface of this cylinder be described by the equation 

x 2 {z-a-Bf 

A 2 + B 2 ' 1 ' 

where A > B are semiaxes and a is the closest separation distance between the lens and 
the plate z = (see Fig. 1). The axis of a cylinder is aligned along the y axis. The upper 
surface of a plate made of the same material coincides with the plane z = 0. The elliptic 
lens under consideration is assumed to be attached to a micromachined oscillator (see Fig. 4 
where a circular cylindrical lens is situated below a plate). In Fig. 1 the plate is 



in Ref. 



placed below a cylindrical lens for convenience in calculations. 

From Eq. ([!]) the explicit equation for the lens surface is given by 



z(x) =a + B-\jB 2 -^x 2 . (2) 
It is assumed that a/B < 1 and the lens is sufficiently thick, so that a/h <C 1 as well. 



Applying the PFA in the general, Derjaguirii 
Fig. 1 in the same way as was done in Ref. 



ormulation 



44j to the configuration of 



for a circular cylinder, one arrives to the 



following Lifshitz-type formula for the Casimir force at temperature T: 

F(a, T) = _ ^ kBTL f qi k ± dk± 

1=0 J ° 

y>3S4 + r&) / dxe- 2n *'W. (3) 



X 

n=l 



Here, fcs is the Boltzmann constant, L is the length of the cylinder which is assumed to be 
infinitely large, k± is the projection of the wave vector on the plane z = 0, qi = (k±+£f /c 2 ) 1 ' 2 , 
£j = 2nkBTl/h with I — 0, 1, 2, . . . are the Matsubara frequencies, and z(x) is defined in 
Eq. ([2j). The primed summation means that the term with I = is multiplied by 1/2. 
The reflection coefficients ttm and tte for the two polarizations of the electromagnetic field 
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(transverse magnetic and transverse electric) are given by 

eiqi - h 



mi + h ' 

qi - h 



r T E = r TE (i£ h k ± ) = 1 \ (4) 

qi + h 

where k t = [k\ + £i£f/c?] and Si = e(i&)- 

Bearing in mind that in accordance with Eq. ([2]) 

x = x (z) = ^ ^B 2 - {a + B - z) 2 , (5) 

D 

one can rearrange Eq. ([3]) to the form 

IksTL v-v f°° 
F(a,T) = V" / qik±dk± 

7T z — ' /n 

00 pa+h 

x&™+4e) / da:(z)e- aB «- 

n=l Ja 

> / qik±dk± (6) 



IknTLA A' 



Z=0 



We next introduce dimensionless integration variables 

£ — a 

v = 2aqi, t = nv (7) 

instead of dimensional k± and z, and dimensionless Matsubara frequencies d = 2a£//c. As 
a result, from Eq. fl6]) we arrive at the expression 

i=0 n=l "'w 
phnv/a i a ^ 

x e - ™ / Bm; e _ *dt. (8) 

Jo v/RT^ 

Here, the reflection coefficients in terms of new variables are given by 



r TM = r TM (iQ,v) 



S l V+^V 2 +(e l -\)(f 

v-^v 2 + ( £l - l)Cf 



r TE = r TE {i( h v) = , ' ==, (9) 

v + ^v 2 + (e, - l)Cf 
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where E\ = e(ic(i/2a). 

Within the application regime of the PFA we are looking for the main contribution to the 
expansion of Eq. in terms of the small parameter a/B < 1. To do so we can restrict our 
consideration to the lowest expansion order in a/B of the integrand with respect to t. We 
can also set the upper integration limit of this integral equal to infinity taking into account 
that h 3> a. This leads to 



na,T) = - k f L 2 " E'Et <1<» 

v 1=0 n=l v 

roo poo -t Jj. 

After calculation of an integral with respect to t, and summation with respect to n one 
obtains 



40Fa 2 V^B hi 
x [U 1/2 (r 2 TM e- v ) + U 1/2 (r 2 TE e- v )} , (11) 



where Li n (z) is the polylogarithm function 45|. This is the Lifshitz-type formula for the 



Casimir force between an elliptic cylinder or cylindrical lens and a plate. For a circular 
cylinder A = B = R, and Eq. (llljl coincides with the result derived in Ref. (36 1. 

In the case of an elliptic cylinder and a plate made of an ideal metal, r\ M = ri| B = 1. 
Then at zero temperature Eq. ffTTl) results in 







00 g— n« 

'n 



£ n=l 



Changing the order of integrations and calculating the integrals one obtains 



^-^^m±h (13) 

v n=l 

After calculating the sum, the result is 



,IM' - ir Lhc A 



= (14) 



For a circular cylinder (A = B = R) this leads to a well known result obtained in Ref. 33]. 



S 



Both Eqs. ( ITT]) and (Tl4"j) are approximate, as they are obtained with the help of the PFA. 



Using the same considerations as in Ref. (36], one can conclude that the relative error of 



these equations is approximately 0.3a/ B. For typical experimental parameters B = 100 /im 
and a = 200 nm the resulting error is equal to 0.06%. 

The Lifshitz-type formula for the gradient of the Casimir force between an elliptic cylinder 
and a plate can be obtained by analogy with Eq. ( TTTj) . For this purpose we differentiate 
Eq. ([3]) with respect to a using Eq. fl2]) and arrive at 



oo 



dF(a,T) Ak B TL^ , 2 



o 



1=0 

x yVr^ + r^) / dxe~ 2nqiz{x) . (15) 

n=l ' ^° 



We then repeat the same transformations as were done in Eqs. (I6l- (fl0l in application to 
Eq. ( TT51) and obtain 



dF{a, T) k B TL A 



da 4v/?a 3 y/2aB 



^ poo 



Z=0 

x [Li_ 1/2 (4 M e-") + Li_ 1/2 (4 E e-*)] . (16) 



For a circular cylinder, Eq. ( TTBl) coincides with the result obtained in Ref. [36|. In the 
case of an ideal metal elliptic cylinder and a plate at zero temperature Eq. (TT6l) leads to 

dF lM {a,0) _ 7n 3 Lhc A 
da' 768a 3 7^1' 

The same equation is obtained by differentiation of Eq. f)14p with respect to a. 

The Lifshitz-type formulas ffTTl) and f[T6"|) allow computations of the Casimir force and 
its gradient in the configuration of an elliptic cylindrical lens and a plate made of real 
materials. In so doing different theoretical approaches can be used, such as the Drude and 
plasma model approaches mentioned in Sec. I. For an Au circular cylinder above a plate 
computations of the relative thermal correction to the Casimir force and its gradient as a 
function of separation using the Drude and plasma model approaches were performed in 
Ref. [36J] within the separation range from 150 nm to 5 /im. It was shown that when the 
Drude model is used the magnitude of the relative thermal corection to the Casimir force 
achieves its maximum value 41.6% at a = 2.55 fim. The maximum magnitude of the relative 
thermal correction to the gradient of the Casimir force 52% occurs at a = 3.6 /im. When 
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the plasma model approach is used, the relative thermal correction to the Casimir force 



increases monotonically from 0.016% at 150 nm to 26.7% at a = 5 fim |36j. In the case of an 
elliptic cylinder the respective results for the relative thermal correction remain the same. 
This allows discrimination between the predictions of different theoretical approaches by 
comparing the computation results with the measurement data. 

By using the PFA, one can also obtain a simple expression for the electric force between 
an elliptic cylinder and a plate. An electric force is used to perform calibrations in the 
measurements of the Casimir force. For a potential difference (V — Vq) between an elliptic 
cylinder and a plate (V is the applied voltage and V is the residual potential), the electric 
force calculated similar to the Casimir force is given by 

where €q is the permittivity of the vacuum. For a circular cylinder A = B = R, this formula 
was obtained 34j from the exact expression for the electric force 46j 



where A = \/h 2 — R 2 and h = R + a. Expanding the right-hand side of Eq. ( |T9l) in powers 
of a small parameter a/R, one obtains 

7reQLyR /rr 2 ( ^ 1 a 17 a 2 



FJa) = ^—(V-Vo) Ml + ; . (20) 

2 v / 2a 3 / 2 V 12 R 480 R 2 J K ' 

Thus, for a = 100 nm and R = 100 /im the error in the electric force due to the use of the 
PFA is equal to only 0.008%. 



IV. AN ASYMMETRIC CYLINDRICAL LENS AND A PLATE 

The above results can be used to calculate the Casimir force between an asymmetric 
cylindrical lens modeled by the two elliptic cylinders with dissimilar semiaxes A±, B\ and 
A 2 , B 2 (see Fig. 2). One half of such a lens of width d\ is produced as a section of an elliptic 
cylinder with semiaxes Ax, B x , and another half of width d 2 as a section of a cylinder with 
semiaxes A 2 , B 2 . In so doing both halves are equal in thickness. Bearing in mind that the 
PFA is an additive method, the Casimir force between each of the halves of an asymmetric 
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cylindrical lens and a plate can be calculated using Eq. ( ITT]) . The Casimir force between the 
entire lens and a plate is then given by 

k B TL 1 / A x A 2 \ 

^--sv^WaKw^wJ (21) 



°° / POO 

xJ2 vV2dv [ L M r ™0 + U 1/2 (r 2 E e 
1=0 



This equation is valid under the conditions ajB\ <C 1, a/B 2 <C 1, and a//i<l. In a similar 
way, the gradient of the Casimir force between an asymmetric elliptic lens shown in Fig. 2 
and a plate is expressed by the equation 

dF(a,T^_ k B TL 1 ( A x , A 2 \ 



da 8v^Fa 3 \^B[ JB 2 

OO 

^ / v 5 / 2 dv[U_ 1/2 (4 M e- v ) + U_ 1/2 (r 2 TE e-n]- 



OO 

X 

1=0 



'0 

Equations (12 lj) and (1221) allow calculation of the Casimir force and its gradient in the config- 
uration of a plate and a cylindrical lens consisting of two parts of dissimilar elliptic cylinders. 

We turn next to the consideration of another asymmetric cylindrical lens which is obtained 
from an elliptic cylinder defined in its proper coordinates (x, z) as a cross section by the 
plane perpendicular to the plane xz and inclined at an angle <p to the axis x [see Fig. 3(a)]. 
We then rotate the resulting lens through an angle <p clockwise around the axis of a cylinder 
in order to make its base parallel to the plate [see Fig. 3(b)]. As before, the thickness of a 
lens is h. 

It is easily seen that the transformation from the coordinates (x, z) to (x, z) shown in 
Fig. 3(b) has the form 

x = Xq + a sin tp + x cos (p — z sin (p, 

z = z — a cos <p + x sin p + z cos <p, (23) 

where (xo, Zq) are the coordinates of the lens point closest to the plate given by 

A 2 sin 09 _, B 2 cos p> , nt . 
^0 = —^, z = ^t. (24) 

Here, we have introduced the notation 



H = H(A, B; <p) = J A 2 sin 2 <p + B 2 cos 2 <p. (25) 
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Now we substitute Eqs. ( 123|) and IH)) into the equation of a lens surface 



written in the proper coordinates and arrive at 

2 A 2 - B 2 . . , A 2 B 2 

x — 2x — — (a — z) sin y? cosy? + 2(a — z) 

+ (a- Z ) 2 A2cOs2 ^f S[n ^ =0. (27) 
H 

This equation describes the surface of an asymmetric cylindrical lens in the coordinates 
(x, z). If the inclination angle is <p = 0, Eq. f[2"T|) simplifies to 

x 2 + 2(a-z)^ + (a- z) 2 ^ = 0, (28) 

and has the solution (jSj) as it must. For a circular cylinder A = B = R and an arbitrary 
angle y?, Eq. ( |2"7|) simplifies to 

x 2 + 2(a - + (a - ;z) 2 = 0, (29) 

leading again to the specific case of Eq. (j3J). 
Equation (|27|) has the following two solutions: 

A 2 - 5 2 

£i,2 = |p — (-2 - a) sin yj cos y? 

± — [2{z - a)H - {z - a) 2 ] l/ \ (30) 

where the upper and lower signs are for x > and x < 0, respectively [see Fig. 3(b)]. 

We next consider the calculation of the thermal Casimir force between an asymmetric 
cylindrical lens and a plate shown in Fig. 3(b). This can be done by using the first equality 
in Eq. ([6]) which we apply separately to the parts of the lens with x < and x > 0: 



u TT °° i f°° °° 
71 1=0 ^° n=l 



X 



pa+h 

dx 2 (z)e- 2nqiz + / dx x (z)e- 2nmz 



.J a+h 

From Eq. (130]) . the differentials dx\£ are given by 



(31) 



(A 2 — B 2 ) sin y? cosy? , 
cteio = —r. dz 

±^ ff "^ + Q 17^ (32) 

# 2 [2(2 -a)F- (z-a) 2 ] 1/2 
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with the same sign convention as formulated above. Substituting Eq. (1321) into Eq. (l3Tj) . one 
obtains 

F y (a, T) = _ 2keTL ^E. f qi kj_dk± (33) 
x^(r TM + r TE) y a [2{z _ a)H _ {z _ a)2] yr 



Introducing the integration variable v from Eq. ([7]) instead of the variable k±, this can be 
rearranged to 



W) = -^§E' /I, 1 * (34) 



47ra 3 if 2 



G 



X y^( r TM + r TE) / 
n=l "' a 



(H-z + a)e- nvz / a dz 



L [2(z-a)H- (z-a) 2 } 1 ' 2 ' 
Now, instead of the variable z, we introduce the variable t defined in Eq. ([7]) and use the 
conditions 

a a , 

- h «l, « 1. (35) 



Then Eq. (jMj) reduces to 



F(a T) = - kBTL ^f' 

x E "7^™ + r Tl)e~ nv / -^dt. (36) 
„_i v n Jo VI 



Performing the integration with respect to t and the summation over n, we arrive at 

F v a,T = --4—-—= — V / v^ 2 dv 

x [Li 1/2 (r 2 M e-) + Li 1/2 (r 2 E e^)]. (37) 

The comparison of this result with Eq. fill I) shows that the dependence of the Casimir force 
on if is contained exclusively in the factor 

/ B\ V2 (A 2 V 3/4 
G(A, B- <p) = I- = I— sin 2 ip + cos 2 , . (38) 

Thus, a similar result is obtained for the gradient of the Casimir force between an asymmetric 
cylindrical lens and a plate 

dF v (a,T) dF(a, T) 

—^—-G(A,B,p)—^, (39) 
13 



where the gradient of the Casimir force between a symmetric elliptic cylindrical lens and a 
plate, dF(a,T)/da, is given by Eq. ( JIBjl . 

From Eq. ( 13 8 p it is seen that the function G and, thus, the Casimir force F v and its 
gradient satisfy the condition 

A g(a,B-(p + ^ = 4=G(B, A; <p). (40) 



'B V 2/ ,/,! 

Specifically, from Eq. f l4"0l) we have 



A G(A,B-*)=* (41) 



1? V 2/ VI 

i.e., the rotation through an angle tp — 7r/2 interchanges the semiaxes of a cylinder, as it 
should. 

In Fig. 4(a) we present the relative Casimir force and its gradient 

M^ = W^ =G(AB ;,) (42) 
F(a,T) oF(a,T)/oa 

as a function of the rotation angle. Different lines are for different values of the ratio of 
semiaxes A/B = 1.1, 1.2, 1.3, and 1.4 increasing from the top to bottom lines. Keeping in 
mind that experimentally it is difficult to ensure exactly ip = 0, we also present in Fig. 4(b) 
the same lines over a narrow interval from <p = to ip = 0.1 rad. From Fig. 4(a) it is seen 
that the relative Casimir force and its gradient decrease monotonically with the increase of (p 
and A/B. According to Fig. 4(b), even the rotation of an elliptic lens through 0.1 rad (5.73°) 
leads to less than 1% deviation of the Casimir force and its gradient from their respective 
values at <p = for any value of A/B considered. A deviation of the Casimir force and its 
gradient from their values at (p — for less than 0.1% is achieved for the rotation angles 
(p < 0.025 rad = 1.43°. This places experimental limitations on an allowed asymmetry of 
the elliptic cylindrical lens used. 

Note that the electric force between an asymmetric cylinder and a plate can be obtained 
from Eq. (ITS]) by the replacement 

* -+1(4= + 4=1 (43) 



'B 2 \^B[ VB 2 / 
or by multiplying the right-hand side of Eq. ( TTSl by the factor G defined in Eq. (1381) de 
pending on the nature of asymmetry. 
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V. MICROMACHINED OSCILLATOR WITH AN ATTACHED CYLINDRICAL 
LENS IN A NONLINEAR REGIME 



In Ref. [36] it was proposed to perform dynamic measurements of the Casimir interaction 
between a plate and a circular cylinder attached to a micromachined oscillator. The proposed 
experiment aims to achieve the same high experimental precision, as in the experiment of 
Refs. 



15 



161 ] for a sphere above a plate, over a wider separation region. For this purpose, the 



same measures, as in Refs. 



15 



161 ] would be undertaken, specifically, to reduce mechanical 



vibrations. At any rate, the effect of vibrations in the position at the proposed measurement 
frequency (a few hundred hertz) is much smaller than the uncertainty in the position due to 
the interferometric technique used. As a result, the impact of vibrations on the determination 
of the gradient of the Casimir force is smaller than the current systematic experimental error, 
and thus, can be neglected. 



As in Refs. 
Ref. 



15 



171 . 1471 . |48| which exploited the configuration of a sphere near a plate, 



discussed measurements of the gradient of the Casimir force in a linear regime (the 
same regime was employed in dynamic measurements by means of an atomic force micro- 



scope 



49 



53|). Here, we find the frequency shift of an oscillator, caused by the Casimir 



force between an elliptic cylinder and a plate, in the nonlinear regime. This allows mea- 
surements down to shorter separation distances where the micromachined oscillator behaves 
nonlinearly. 

In the dynamic regime the separation distance between an elliptic cylinder attached to a 
micromachined oscillator and a plate is varied with time harmonically 



a(t) = a + A z cosu r t. 



(44) 



Here, u r is the resonant frequency of the oscillator under the influence of the Casimir force 
acting between a cylinder and a plate. The amplitude of oscillations A z should be sufficiently 
small in comparison with the separation a. In the presence of the Casimir force, F(a,T), 



the frequency u r is different from the natura 
an oscillator problem was considered in Refs. 



angu 



54 



ar frequency of the oscillator uq. Such 



55| perturbatively and in Ref. |56| exactly. 



The exact expression for the shift of the second power of the natural frequency of an oscillator 
produced by the Casimir force is given by 



Cu> r 

ixA z 



dt cos(oj r t)F(a + A z coscovt, T). 
15 



(45) 



Here, C is a constant depending on specific parameters of the setup used. Thus, for a 
micromachined oscillator C = b 2 /I, where b and I are the lever arm and the moment 



of inertia. Note that in Ref. 56], where the Bose-Einstein condensate was considered as 
a second body, the Casimir-Polder force between individual atoms and a plate was also 
averaged over the condensate cloud. 

In the case of an elliptic cylinder interacting with the plate the force F is given by 
Eq. (ITT)) . Representing the polilogarithm functions in Eq. (ITT)) as power series, and replacing 
the integration variable v with k± = \/v 2 — Qj (2a), we rearrange this equation to the form 

k B TL A ™>™ 1 [°l 

v v ^ i=o n=l v Jl) 

x (4 n M + 4 n E)e" 2a9in . (46) 

Substituting Eq. ( H6j) into Eq. ( 145)) and introducing the new integration variable 9 = u r t, 
we arrive at 

a 2 c k B TL a in., r- 

/■27T 

x (4m + 4 n E ) e - 2 ^ n / d# cos 6e- 2nAzqi cose . (47) 

Integrating by parts, the integral with respect to 6 can be reduced to [57] 

sin^e ±2cos ^= | (48) 

where /fc(-z) is the Bessel function of an imaginary argument. The frequency shift (ITT)) is 
then expressed as 



a 2 2C k B TL A 1 
w r - w = — p= = > > — = / k ± dkj 



x Voi(r^ M + r 2 T n E )e- 2a «' n h(2A zqi n). (49) 



Returning to the variable v, one finally obtains 

C k B TL A 



OO ^ OO j 

777 



A 2 2v7Fa 2 v^TB ^ ^ v 7 ^ 



x J v 3/2 dv(r^ n M + r^)e™h [^nvj . (50) 

This is the general expession for the shift of the second power of the oscillator frequency due 
to the Casimir force between an elliptic cylindrical lens and a plate obtained using the PFA. 

16 



Equation (!50|) takes into account nonlinearity of the oscillator. By making the replacement 
(|43p it can be generalized to the case of an asymmetric cylindrical lens consisting of the 
parts of two dissimilar elliptic cylinders. Multiplying the right-hand side of Eq. ( I50p by the 
factor G(A, B; (p) defined in Eq. (1551) one obtains the generalization of Eq. ([5"U]) for the case 
of elliptic cylindrical lens rotated through an angle tp (see Sec. IV). 

It is easily seen that in the linear approximation Eq . ([50 ) leads to familiar expressions 

551 ] . In fact the linear regime of 



15. 



16 



47 



48 



54 



commonly used in the literature 

the oscillator, and the first nonlinear corrections to it considered in Ref. 
by using the following representation for the Bessel function 



551 ]. are obtained 



57j 



A(*) = £ 



(I) 



2fc+l 



fc=0 



k\(k + l)\ 



2 16 V 1 



(51) 



Substituting the first term on the right-hand side of this equation into Eq. (150|) and per- 
forming the summation over n, one finds 

k B TL A 



-c 



4v^a 3 ^/2a~B 



J2 / V^d 

1=0 



X 



[Li_ 1/2 (r| M e ")|LL 1/2 (4e ")] . 



(52) 



We can now use t 
arrive at 



-1U w u se I 



ne 



16 



ifshitz-type formula (flBT) for the gradient of the Casimir force, and 



47|,|48 



54 



-c 



dF(a,T) 
da 



(53) 



Bearing in mind that and, thus, co r + u Q ~ 2u , Eq. (|53|) is often presented in the 

form 

C dF(a,T)' 



0J r 



U) 



2ul 



da 



(54) 



In the linear regime, Eqs. ( 153]) and (|54j) allow calculation of the frequency shift of the 
oscillator due to the Casimir force between a plate and an elliptic cylinder. However, beyond 
the linear regime, the frequency shift should be calculated using Eq. ( 150]) . This allows reliable 
comparison between the experimental data and theoretical results for the Casimir force. 



VI. CONCLUSIONS AND DISCUSSION 



In this paper we have investigated the Casimir force acting between a plate and a micro- 
fabricated elliptic cylindrical lens made of real materials. This problem is of topical interest 

17 



from the experimental point of view. The application of available technologies discussed by 
us in Sec. II leads to a fabrication of elliptic cylindrical lenses on top of a micromachined 
oscillator, rather than just the circular cylindrical lenses considered previously in the litera- 
ture. We have obtained the Lifshitz-type formulas for the Casimir force and for the gradient 
of the Casimir force in the configuration of an elliptic cylinder and a plate. In the framework 
of the PFA (which is applicable at separations much less than the smaller semiaxis of an 
elliptic cylinder), the results for an elliptic cylinder are obtained from the respective results 
for a circular cylinder by replacing the cylinder radius R with A 2 /B, where A and B are 
the semiaxes of an elliptic cylinder (see Sec. III). 

Bearing in mind that nanotechnological fabrication procedures may lead to cylinders 
with deviations from perfect elliptic shape, we considered two types of such deviations. In 
Sec. IV we obtained the Lifshitz-type formulas for the Casimir force and for its gradient 
in the configuration of a plate near asymmetric elliptic cylindrical lenses. Specifically, the 
constraints on an admissible angle of rotation of an elliptic cylindrical lens about the cylinder 
axis were found, allowing sufficiently small deviations from the values of the Casimir force 
and its gradient computed for the case of zero rotation angle. The respective results for both 
perfect and asymmetric elliptic cylindrical lenses were also obtained for the electrostatic 
force in plate-lens configuration used for calibration purposes in experiments on measuring 
the Casimir force. Note also that corrections to the Casimir force and its gradient due to 
nonparallelity of a plate and an elliptic cylinder are approximately the same as in the case 
of circular cylinder considered in Ref. 36]. 

For the needs of several proposed experiments on measuring the Casimir force in a 
cylinder-plate geometry, we have considered an oscillator with an attached elliptic cylin- 
drical lens interacting with the plane plate both made of real materials. For dynamic mea- 
surements, when the separation distance between a lens and a plate is varied harmonically, 
we have found the frequency shift of an oscillator due to the Casimir force in the nonlinear 
regime (Sec. V). The resulting equations can be used at short separations between a lens 
and a plate where the commonly used linear equations are not applicable. At the same time, 
it is shown that in the linear approximation our result yields to the known expression. 

To conclude, the proposed experiment on measuring the Casimir force between a micro- 
fabricated elliptic cylindrical lens on the top of a micromachined oscillator and a plate is 
of much current interest and can shed additional light on the problem of thermal Casimir 



18 



force. 



Acknowledgments 

R.S.D. acknowledges NSF support through Grant No. PHY-0701236 and LANL support 
through contract No. 49423-001-07. D.L. and R.S.D. acknowledge support from DARPA 
grant No. 09-Y557. E.F. was supported in part by DOE under Grant No. DE-76ER071428. 
G.L.K. and V.M.M. are grateful to the Department of Physics, Purdue University for finan- 
cial support. G.L.K. was also partially supported by the Grant of the Russian Ministry of 
Education P-184. 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. B 51, 793 (1948). 

[2] M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Advances in the 

Casimir Effect (Oxford University Press, Oxford, 2009). 
[3] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Rev. Mod. Phys. 81, 1827 

(2009). 

[4] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Int. J. Mod. Phys. B 25, 171 
(2011). 

[5] G. L. Klimchitskaya, J. Phys.: Conf. Series 161, 012002 (2009). 
[6] M. Bostrom and B. E. Sernelius, Phys. Rev. Lett. 84, 4757 (2000). 
[7] J. S. H0ye, I. Brevik, J. B. Aarseth, and K. A. Milton, Phys. Rev. E 67, 056116 (2003). 
[8] I. Brevik, J. B. Aarseth, J. S. H0ye, and K. A. Milton, Phys. Rev. E 71, 056101 (2005). 
[9] M. Bordag, B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. Lett. 85, 
503 (2000). 

[10] C. Genet, A. Lambrecht, and S. Reynaud, Phys. Rev. A 62, 012110 (2000). 
[11] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, J. Phys. A: Math. Theor. 40, 
339 (2007). 

[12] V. B. Bezerra, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. A 65, 052113 
(2002); ibid 66, 062112 (2002). 



19 



[13] V. B. Bezerra, G. L. Klimchitskaya, V. M. Mostepanenko, and C. Romero, Phys. Rev. A 69, 
022119 (2004). 

[14] M. Bordag and I. Pirozhenko, Phys. Rev. D 81, 085023 (2010). 

[15] R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Phys. Rev. D 75, 077101 (2007). 

[16] R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Eur. Phys. J. C 51, 963 (2007). 

[17] R. S. Decca, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, D. Lopez, and V. M. Mostepa- 
nenko, Phys. Rev. D 68, 116003 (2003). 

[18] B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. A 82, 032513 (2010). 

[19] F. Intravaia and C. Henkel, Phys. Rev. Lett. 103, 130405 (2009). 

[20] G. Bimonte, Phys. Rev. A 79, 042107 (2009). 

[21] B. E. Sernelius, Phys. Rev. A 80, 043828 (2009). 

[22] P. R. Buenzli and Ph. A. Martin, Phys. Rev. E 77, 011114 (2008). 

[23] V. M. Mostepanenko and G. L. Klimchitskaya, Int. J. Mod. Phys, A 25, 2302 (2010). 

[24] T. Emig, R. L. Jaffe, M. Kardar, and A. Scardicchio, Phys. Rev. Lett. 96, 080403 (2006). 

[25] O. Kenneth and I. Klich, Phys. Rev. Lett. 97, 160401 (2006). 

[26] M. Bordag, Phys. Rev. D 73, 125018 (2006). 

[27] T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. Rev. Lett. 99, 170403 (2007). 

[28] A. Lambrecht and V. N. Marachevsky, Phys. Rev. Lett. 101, 160403 (2008). 

[29] H.-C. Chiu, G. L. Klimchitskaya, V. N. Marachevsky, V. M. Mostepanenko, and U. Mohideen, 

Phys. Rev. B 81, 115417 (2010). 
[30] A. O. Sushkov, W. J. Kim, D. A. R. Dalvit, and S. K. Lamoreaux, Nature Phys. 7, 230 (2011). 
[31] V. B. Bezerra, G. L. Klimchitskaya, U. Mohideen, V. M. Mostepanenko, and C. Romero, Phys. 

Rev. B 83, 075417 (2011). 
[32] M. Masuda and M. Sasaki, Phys. Rev. Lett. 102, 171101 (2009). 

[33] D. A. R. Dalvit, F. C. Lombardo, F. D. Mazzitelli, and R. Onofrio, Europhys. Lett. 67, 517 
(2004). 

[34] M. Brown-Hayes, D. A. R. Dalvit, F. D. Mazzitelli, W. J. Kim, and R. Onofrio, Phys. Rev. 
A 72, 052102 (2005). 

[35] Q. Wei, D. A. R. Dalvit, F. C. Lombardo, F. D. Mazzitelli, and R. Onofrio, Phys. Rev. A 81, 

20 



052115 (2010). 

[36] R. S. Decca, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, D. Lopez, and V. M. Mostepa- 

nenko, Phys. Rev. A 82, 052515 (2010). 
[37] L. P. Teo, Phys. Rev. D 84, 025022 (2011). 

[38] V. B. Bezerra, E. R. Bezerra de Mello, G. L. Klimchitskaya, V. M. Mostepanenko, and 

A. A. Saharian, Eur. Phys. J. C 71, 1614 (2011). 
[39] S. Reyntjens and R. Puers, J. Micromech. Microeng. 11, 287 (2001). 

[40] S. Matsui, Three-dimensional nanostructure fabrication by focused ion beam chemical vapor 
deposition, in Springer Handbook of Nanotechnology, ed. B. Bhushan (Springer, Heidelberg, 
2010), p.211. 

[41] Y. Cheng, H. L. Tsai, K. Sugioka and K. Midorikawa, Appl. Phys. A 85, 11 (2006). 

[42] S. Campbell, The Science and Engineering of Microelectronic Fabrication, 2nd edn. (Oxford 

University Press, 2001). 
[43] S. Franssila, Introduction to Micro fabrication (John Wiley &; Sons, London, 2004). 
[44] B. V. Derjaguin, Kolloid. Z. 69, 155 (1934). 

[45] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, Integrals and Series, Vol. 2 (Gordon 
and Breach, New York, 1986). 

[46] W. R. Smythe, Static and Dynamic Elictricity (McGraw-Hill, New York, 1968). 

[47] G. L. Klimchitskaya, R. S. Decca, E. Fischbach, D. E. Krause, D. Lopez, and V. M. Mostepa- 
nenko, Int. J. Mod. Phys. A 20, 2205 (2005). 

[48] R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Ann. Phys. (NY.) 318, 37 (2005). 

[49] G. Jourdan, A. Lambrecht, F. Comin, and J. Chevrier, Europhys. Lett. 85, 31001 (2009). 

[50] S. de Man, K. Heeck, R. J. Wijngaarden, and D. Iannuzzi, Phys. Rev. Lett. 103, 040402 
(2009). 

[51] S. de Man, K. Heeck, and D. Iannuzzi, Phys. Rev. A 82, 062512 (2010). 

[52] G. Torricelli, P. J. van Zwol, O. Shpak, C. Binns, G. Palasantzas, B. J. Kooi, V. B. Svetovoy, 

and M. Wuttig, Phys. Rev. A 82, 010101 (R) (2010). 
[53] G. Torricelli, I. Pirozhenko, S. Thornton, A. Lambrecht, and C. Binns, Europhys. Lett. 93, 

51001 (2011). 

[54] H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Science, 291, 1941 

21 



(2001). 

[55] H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Phys. Rev. Lett. 
87, 211801 (2001). 

[56] M. Antezza, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 70, 053619 (2004). 
[57] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products (Academic Press, 
New York, 1980). 



22 




FIG. 1: Elliptic cylindrical lens of thickness h and width 2d obtained from an elliptic cylinder with 
semiaxes A and B above a plate. The figure is not to scale. 
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FIG. 2: An asymmetric elliptic cylindrical lens of thickness h and width d\ + d,2 obtained from two 
elliptic cylinders with semiaxes A\, B\ and A2, B2 above a plate. The figure is not to scale. 
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FIG. 3: (a) An asymmetric elliptic cylindrical lens of thickness h obtained from an elliptic cylinder 
with semiaxes A and B. (b) The same asymmetric elliptic cylindrical lens spaced at the closest 
separation a above a plate. The figure is not to scale. 
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FIG. 4: (Color online) The ratio of the Casimir forces between an asymmetric elliptic cylindrical 
lens and a plate for the angle of rotation equal to <p and to zero as a function of (p. For different 
lines the ratio of cylinder semiaxes A/B = 1.1, 1.2, 1.3, and 1.4 increasing from the top to bottom 
lines. The interval of the angles of rotation varies (a) from ip = to (p = tv/2 and (b) from <p = 
to <p = 0.1 rad. 
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